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ABSTRACT 

In this Letter, we analyse the predicted physical properties of massive galaxies, in the 
framework of recent semi-analytic models of galaxy formation. All models considered 
account for winds driven by supernovae explosions and suppression of gas condensation 
at the centre of relatively massive haloes by active galactic nuclei (AGN). We show 
that, while these models successfully reproduce the old stellar populations observed 
for massive galaxies, they fail in reproducing their observed chemical abundances. 
This problem is alleviated but still present if AGN feedback is completely switched 
off. Moreover, in this case, model predictions fail in accounting for the old stellar 
ages of massive galaxies. We argue that the difficulty of semi-analytical models in 
simultaneously reproducing the observed ages and metallicities of massive galaxies, 
signals a fundamental problem with the schemes that are currently adopted to model 
star formation, feedback, and related recycling of gas and metals. 
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1 INTRODUCTION 

In the current standard picture for structure formation, 
galaxies originate from gas condensation within the poten- 
tial wells of hierarchically formed dark matter haloes. The 
mass distribution of these haloes is well described by a 
power-law, differing significantly from the typical shape of 
the observed galaxy mass function. 

Stellar feedback i s believed to play a c rucial role at 
the l ow mass regime (|White fc Frenkl Il99ll : iBenson et alj 
2003). Observations suggest that galactic-scale outflows are 
ubiquitous in starburst galaxies at all cosmic epochs, and 
that the outflowing material is multiphase (containing cold, 
warm, and hot gas, dust and magn eti zed re lativistic plasma 
- lHeckmanl e.g. l2002l : IWeiner et al.l e.g. l2009l ). Unfortunately, 
available observational measurements refer to material that 
is still relatively deep within the gravitational potential of 
the halo. So the estimated outflow rates are difficult to 
translate into rates at which mass, metals, and energy es- 
cape from galaxies and are eventually transported into the 
inter-galactic medium. The fate of the outflowing material 
will depend critically on a number of unknowns, as well as 
on its multiphase nature. Given the uncertainties, it is not 
clear how appropriate the different prescriptions adopted for 
treating galactic winds in galaxy formation models are. 

At the massive end, feedback from active galactic nu- 
clei (AGN) is believed to play a key role. X-ray observa- 
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tions show that AGN feedback should be responsible for the 
thermal structure of "cool cores" ; only a modest amount of 
gas cools and form stars at the centre of galaxy clusters, 
despite the short cooling timescales inferred from the X-ray 
emission. Detailed recent studies have confirmed that bright- 
est cluster galaxies (BCGs) are more likely to host a radio- 
loud AGN than other galaxies of similar stellar mass, and 
have shown that the ensemble averaged power from radio 
galaxies is suffic i ent to offset the mean level of cooling (e.g., 
iBest et al.ll2007l : iMcNamara fc Nulsenll2007l ). Heating from 
a central AGN has become a crucial ingredient for galaxy 
formation models in order to reproduce the observed expo- 
nential cut-off at the bright end of the galaxy luminosity 
function, and the old stellar populations observed for mas - 
sive galaxies (e.g. ICroton et ai1l2006l ; iDe Lucia et alj|2006h . 
The details of the energy injection by the central engine, and 
coupling with the surrounding hot gas are, however, still un- 
clear. So the prescriptions adopted to model this process are 
very schematic, and often not well grounded in observations. 



2 THE GALAXY FORMATION MODELS 

In this study, we take advantage of publicly available galaxy 
catalogue^, constructed applying semi-analytic methods 
to high-resolution cosmological simulations. In particular, 



1 These are available at the following 
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Bower et al.1 l|2006l. 
20071 . DLB07 ), and 



we use catalogues from the models by 
BOW06 he r eafter ) , |Pe Lucia fc Blaizotl 

iGuo et al.1 J201ll. GUOll ) applied to the Millennium Simu- 
lation i|Springel et al.ll2005r ). 

The first two models have been developed indepen- 
dently and differ in a number of details related to the merger 
tree construction and the prescriptions adopted for various 
physical processes considered. The GUOll model is based on 
the DLB07 one, but it adopts a different (significantly more 
efficient) supernovae (SN) feedback and an improved treat- 
ment for the evolution of satellite galaxies. All three models 
include prescriptions for SN driven winds, follow the growth 
of super-massive black holes, and include a phenomenolog- 
ical description of AGN feedback. All three models adopt 
an instantaneous recycling approximation, and assume an 
efficient mixing of metals within the cold gaseous phase. 

Using the public catalogues, we have selected all galax- 
ies in a sub-cube of the Millennium Simulation of 100 Mpc 
on a side. In addition, we have selected all central galaxies of 
haloes more massive tharQ Afaoo = 5 x 10 14 M . We focus, 
in particular, on the stellar ages and metallicities of model 
galaxies. Unfortunately, the public catalogues do not provide 
consistent information: ages and metallicities are weighted 
by mass in the DLB07 and GUOll models, and by V-band 
luminosity for the BOW06 model. This does not affect, how- 
ever, our conclusions. In the following, we will assume for 
the solar metallicity: Zq = 0.02. 



3 THE STELLAR POPULATIONS OF MODEL 
GALAXIES 

The gray regions in Figure [1] show the age-stellar mass re- 
lation for all model galaxies in the sub-volume of the Mil- 
lennium simulation considered. Solid and dashed lines show 
the median and the 16th and 84th percentiles of the dis- 
tributions, while filled symbols show the location of central 
galaxies of haloes with M200 > 5 x 10 14 Mq. As mentioned 
above, the stellar ages from the DLB07 and GUOll models 
are weighted by mass, while those from the BOW06 model 
are weighted by V-band luminosity. The latter are more 
directly comparable to the observational measurements by 
iGallazzi et all (2005), shown as purple lines. These are ob- 
tained adopting a Bayesian statistical approach to derive full 
likelihood distributions for ages and metallicities by compar- 
ing the observed spectra to a library of model star formation 
histories. The comparison is based on the strengths of spec- 
tral absorption features that depend weakly on the [a/Fe] 
ratio. The distributions obtained for model galaxies are not 
convolved by observational uncertainties, that are relatively 
large particularly for older ages. 

Model galaxies tend to be on average older than ob- 
servational estimates, particularly for low-mass galaxies. 
This problem is in part due to the wel l known excess of 
faint and passive mod el satellites (e.g. IWang et al.l I2007J; 
IWeinmann et al1l2010l . and references therein), ft appears 
less severe when considering luminosity-weighted ages (see 
right panel), but models still miss a population of very young 



2 M200 is the mass contained within the radius encompassing an 
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low-mass galaxies. At the massive end, model galaxies are 
dominated by very old stellar populations so that differences 
between mass-weighted and luminosity-weighted ages are 
small. In this mass range, model predictions are close to the 
upper limits of the observational estimates. Central galaxies 
of massive clusters have very old ages: in the DLB07 model, 
the median age is 11.8 Gyr with a scatter of only about 
0.3 Gyr. For the same galaxies, the GUOll model predicts 
somewhat younger mass-weighted ages (the median age is 
11.4 Gyr) with a comparable scatter. This is likely due to 
the fact that, in this model, the hot reservoir associated 
with satellite galaxies is not stripped instantaneously. Cool- 
ing is allowed on these galaxies so that they are on average 
more gas-rich than in a model with instantaneous stripping 
of hot gas. Their accretion onto the central galaxies trigger 
star formation episodes that slightly rejuvenate the stellar 
population of the remnant galaxies. For these galaxies, the 
median luminosity weighted age from the BOW06 model is 
11.2 Gyr, with a ~ 0.4 Gyr scatter. 

Figure [2] shows model predictions and observational es- 
timates for the stellar metallicity-mass relation, with lines 
and symbols having the same meaning as in Figure [1] All 
models predict a relatively tight metallicity-mass relation 
with a steep slope and a pronounced turn-over at the most 
massive end. Central galaxies of massive clusters have stel- 
lar metallicity of ~ 0.57 Zq in the DLB07 model, with a 
~ 0.04 scatter. Because of the more efficient supernovae 
feedback, the GUOll model predicts a steeper metallicity- 
mass relation. For the central massive galaxies considered, 
this model predicts a median metallicity of ~ 0.74 Zq with 
a ~ 0.08 scatter. The luminosity weighted metallicities from 
the BOW06 model are offset low with respect to predictions 
from other models, as expected. For the central galaxies of 
massive clusters, this model predicts a median metallicity of 
~ 0.43 Z©, with a ~ 0.10 scatter . Purple lines show mea- 
surements by IGallazzi et all |2005). The typical uncertainty 
on stellar metallicity is ~ 0.12 dex, but there is a signifi- 
cant tail of galaxies with uncertainties of up to ~ 0.25 dex. 
Low-metallicity galaxies tend to be associated with larger 
uncertainties because of the weaker absorption lines. 



4 THE STELLAR POPULATIONS OF BCGS 

The models used in this study successfully reproduce the 
old stellar populations observed for massive galaxies. If any, 
model galaxies are too old with respect to data. All three 
models considered, however, fail in reproducing the observed 
chemical abundances of massive galaxies: at the massive end, 
the predicted metallicity mass relation turns-over and is off- 
set low with respect to the data. The problem appears to be 
more severe when considering the most massive galaxies at 
the centre of relatively massive haloes, but is not limited to 
them. This is due to the fact that many of the most massive 
satellit es have been accreted recently from relatively massive 
haloes (|De Lucia et al.ll2012l ). Therefore, the same physical 
processes acting on today's BCGs have played an important 
role during the lifetime of these galaxies. 

From the observational viewpoint, while stellar popula- 
tions in early type galaxies have been extensively studied, 
very little is kn own about the stellar popula tions of BCGs. 
In a recent work. lvon der Linden et ahl (|2007h studied a sam- 
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Figure 1. Age-stellar mass relation for the three models considered. Grey maps show the distribution of model galaxies in a sub- volume 
of the Millennium Simulation of size 100 Mpc on a side. Red solid lines show the median, while dashed red lines mark the 16th and 84th 
percentiles of the distributions. Filled circles show the location of the central galaxies of haloes with M200 > 5 X 10 14 Mq. The DLB07 
and the GUOll model ages are weighted by mass, while the BOW06 model ages are weigh t ed by V-band luminosity. Purple lines show 
the conditional distribution of stellar ages as a function of stellar mass from lGallazzi et al.l [j2005T ). 




pie of ~ 600 BCGs and contrasted their stellar populations 
to those of a 'control-sample' of non-BCGs matched in stel- 
lar mass, redshift and colour. They found that the ages and 
the metallicities of BCGs are comparable to those typical 
of massiv e ellipticals, wh i le the a-enhancements are higher 
in BCGs. lLoubser et al.l (|2009h analysed the stellar popu- 
lations of ~ 50 nearby BCGs using long-slit spectroscopy. 
They found that most BCGs are very old, but some (gener- 
ally associated with cooling flows clusters) show signatures 
of recent star formation. The metallicity and Q-enhancement 
distributions measured by Loubser et al. peak at values 
higher than those obtained for normal giant ellipticals. It is 
worth noting that these measurements apply to the central 
regions of BCGs: the inner 3 arcsec in the study by von der 
Linden et al., and 1/8 of the half-light radius in Loubser et 



al. However, the mean stellar population gradients of BCGs 
(null for a-enhancements, almost null for age, and negativ^3 
for metallicity) are co nsistent with those of normal mas - 
sive elliptical galaxies (jLoubser fc Sanchez-Blazquezl 120121 ) ■ 
Therefore, any attempt to correct for aperture effects would 
change the normalization but not the shape of the model 
mass-metallicity relation at the massive end, making the dis- 
crepancy shown in Fig. [2] robust and severe. 

In models where cooling flows are suppressed at late 
times by AGN activity (like those considered here) , the stars 
of massive central galaxies form early, in low-mass progen- 
itors whose high star formation rates are fuelled by rapid 



3 Around -0.3 dex per decade of variation in the radius. 
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Figure 3. Distribution of stellar masses (left), mass- weighted ages (middle), and stellar metallicities (right) for central galaxies of haloes 
with mass M200 > 5 X 10 14 Mq. Black unfilled histograms correspond to the fiducial model described in DLB07 , while red hatched ones 
correspond to the same model but with AGN feedb ack switched off . Purple cross hatched histograms correspond to a model with AGN 
feedback but the SN feedback scheme used in lDe Lucia et aL I d2004) . 



cooling. Galaxies are then assembled relatively late, through 
the accretion of many smaller satellites, driven by the merg- 
ing history of the parent halo. The metallicity-mass relation 
does not evolve significantly in the models, so that these late 
dry mergers tend to decrease their total stellar metallicity. 
Dry minor mergers would likely deposit metal poor stars in 
the outer regions of the remnant, so they would not affect 
significantly the central stellar metallicity of BCGs. For the 
DLB07 model, we find that the stellar metallicity of the main 
progenitor of the BCGs at z ~ 1 is on average only about 
3 per cent lower than its final value, so that even removing 
all stars accreted after z ~ 1 does not improve the disagree- 
ment with observational data. Therefore, in a scenario where 
the mass growth of the most massive galaxies is dominated 
by minor dry mergers, it appears difficult to reproduce si- 
multaneously the observed mass-metallicity relation and the 
chemical abundances of the most massive galaxies. 

It is interesting to analyse what happens when AGN 
feedback is switched off. To this aim, we have re-run the 
DLB07 model on all merger trees of haloes more massive 
than 5 x 10 14 Mq. Fig. [3] shows the distribution of stellar 
masses, mass-weighted ages, and stellar metallicities for all 
central galaxies of these haloes. Black unfilled histograms 
show results from the DLB07 model, while red hatched ones 
correspond to the same model but with no AGN feedback. 
As expected, galaxies become significantly more massive (by 
a factor ~ 6) in a model with no AGN feedback. Gas cool- 
ing and, in minor part, accretion of other gas-rich satellites 
provide fresh material for late star formation, which makes 
the stellar population of these galaxies significantly younger 
(by ~ 5 Gyr). Since the material out of which new stars 
form is enriched by past generations of stars, the mean stel- 
lar metallicity increases going from ~ 0.6 Zq to ~ 0.7 Zq. 
The increase is not large because the metallicity of the cool- 
ing gas is relatively low (~ 0.1 — 0.2 Zq solar, in agreement 
with observational measurements), and it is further diluted 
by metals being deposited t hrough accreted galaxies . Inter- 
estingly, a recent study by iMcCarthv et all (2010) shows 



that hydrodynamical simulations including AGN feedback 
provide central group galaxies that have too low metallicity 
with respect to observational data. They find, however, that 
when AGN feedback is switched off, the metal content of 
central galaxies increases dramatically, becoming even larger 
than observed. 



5 DISCUSSION AND CONCLUSIONS 

All recent models of galaxy formation combine a relatively 
strong SN feedback with a model for AGN feedback to sup- 
press cooling at the centre of relatively massive haloes. The 
former is needed to bring the faint-end slope of the galaxy 
luminosity /mass function in agreement with the relatively 
shallow value measured. The latter affects significantly the 
number densities and stellar ages of the most massive galax- 
ies. Our observational and physical understanding of both 
processes is rather poor so that both are described in galaxy 
formation models using quite schematic prescriptions. 

We have shown that models fail dramatically in repro- 
ducing the observed stellar populations of galaxies. At the 
low mass end, they tend to over-predict the fraction of pas- 
sive galaxies, and lack a population of very young galax- 
ies that is observed. This is a well known problem, that 
plagues all recently published models as w ell as hydrody- 
namical simulations (|Weinmann et ai] l2012). In this letter, 
we have shown that significant discrepancies are found also 
for the most massive galaxies, whose model stellar popula- 
tions are very old but relatively metal poor. 

As mentioned above, the modelling adopted for AGN 
feedback is quite simple and neglects some relevant aspects 
of AGN activity (e.g. a finite duty-cycle, a better modelling 
of t he interaction between the cooling gas and radio jets - see 
also iFontanot et alj|20lj ). Thus, there is room for residual 
star formation at late times in BCGs, which would increase 
the metallicity of their stellar population while decreasing 
their ages. However, we have shown that the stellar metal- 
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licities of the most massive galaxies would be too low with 
respect to data even in the case AGN feedback is switched 
off. An improved treatment of the satellite evolution that al- 
lows residual star formation in satellite galaxies would help 
increasing the stellar metallicity of massive galaxies. How- 
ever, predictions from the GUOll model (that uses such a 
scheme) are still off the observational measurements. 

Uncertainties in the stellar yields can affect model stel- 
lar metallicities. However, the DLB07 and GUOll models 
assume already a relatively large yield (y = 0.03) while the 
BOW06 model assumes y = 0.02 which is closer but still 
somewhat larger than a 'standard' value. In addition, one 
should consider that increasing the yield would increase not 
only the total amount of metals in the stars (and in the gas), 
but also the overall luminosities of model galaxies (because 
of the strong dependence of cooling rates on metallicity) . A 
variable Initial Mass Function (IMF) could also affect the 
total metal content of the most massive galaxies if these 
formed earlier than their lower mass counterparts, from ma- 
terial with different physic al properties. For example, hydro- 
dynamical simulations bv lSmith fc Sigurdssonl (|2007l ) show 
that above a critical metallicity of about 10~ 3 Zq clouds 
can fragment to form low-mass stars, while for gas of lower 
metallicities stars form following a more top-heavy IMF. 
This critical metallicity value, however, is well below that of 
observed galaxies. Finally, our chemical enrichment scheme 
is rather crude, neglecting the mass dependent lifetimes of 
stars, the influence of inefficient mixing and metal loading. 
All these processes can affect the distribution of metals in 
different baryonic components. To illustrate the importance 
of SN feedback, we show in Figure [3] predictions from a 
model with AGN feedback but an alternative SN feedback 
(adopted in lDe Lucia. Kauffmann~fc Whitdl2004 based on 
energy conservation arguments - purple cross hatched his- 
tograms). As explained in DLB07 , this model results in less 
efficient outflows and therefore longer star formation histo- 
ries and higher stellar metallicities. However, it also over- 
predicts the overall luminosities of model galaxies. This ex- 
ample highlights the sensitivity of metallicity on the details 
of SN feedback models. Therefore, requiring a model to pre- 
dict at the same time the correct number densities, ages 
and metallicities of galaxies provides strong constraints on 
the processes responsible for the gas and metal recycling, 
and on the time scales involved. 

In future work, we will investigate how results change 
when more detailed chemical enrichment and recycling 
schemes are adopted. Some fine-tuning might be needed in 
order to avoid making the most massive galaxies too young. 
Therefore, it appears also very important to obtain better 
stellar population estimates for BCGs, and to understand 
how observational estimates compare with the mass and/or 
luminosity weighted values from models. 
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